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Abstract 

The WITCH experiment aims to study a possible admixture of a scalar or tensor type in- 
teraction in /3-decay by determining the P — u angular correlation from the shape of the recoil 
energy spectrum. The installation period was completed and intensive commissioning of the 
set-up was performed already. The lay-out of the WITCH set-up and results of commissioning 
tests performed until now are described here, showing that the full set-up up to the spec- 
trometer is now operational, although several efficiencies are still to be improved. Due to its 
feature of being able to measure the energy spectrum for recoil ions, the WITCH experiment 
also opens possibilities for other observables. 

Key words: weak interaction, beta decay, scalar interaction. Penning trap, retardation spec- 
trometer 



1 Introduction 

Despite the fact that the /3-decay process was discovered already at the end of the 19th century, 
our understanding of weak interactions has developed only gradually. The Standard Model of the 
electroweak interaction is very successful in describing existing experimental data both qualita- 
tively and quantitatively, but a number of parameters have to be determined experimentally and 
several important properties of the interaction are not well understood. The most general interac- 
tion Hamiltonian for nuclear /3-decay which includes all possible interaction types consistent with 
Lorentz-invariance is given by P, Q 

* Corresponding author: [Valentin. Kozlov@fys.kuleuven.ac.be| 
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The coefBcients Ci and C', i £ {S,V,T,A,P} are the coupling constants for Scalar (S), Vector 
(V), Tensor (T), Axial- Vector (A) and Pseudoscalar (P) contributions. In the Standard Model of 
the weak interaction only V and A interactions are present at a fundamental level, which leads to 
the well-known V — A structure of the weak interaction. However, this assumption is based on 
experimental results only and the presence of scalar and tensor types of weak interaction is today 
ruled out only to the level of about 8% of the V- and A- interactions 01 ■ 

A possible admixture of a scalar or tensor type weak interaction in /^-decay can be studied by 
determining the (3 — v angular correlation. This (3 — v angular correlation for unpolarized nuclei 
can be written as 0|: 
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where 9p„ is the angle between the (3 particle and the neutrino, E, vp/c and m are the total energy, 
the velocity relative to the speed of Hght and the rest mass of the /3 particle, F = ^1 — {aZy 
with a the fine-structure constant and Z the nuclear charge of the daughter nucleus, 6 is the 
Fierz interference term which has experimentally been shown to be small (e.g. < 0.0044 at 
90% C.L. 0) and can as a first approximation be assumed to be zero, and a is the (3 — v angular 
correlation coefficient. Since the S- and V-interactions lead to Fermi transitions (F) and the A- 
and T- interactions to Gamow- Teller transitions (GT), the (3 — v angular correlation coefficient a 
can be approximated as (assuming maximal parity violation for the V- and A-interactions and real 
couplings) 
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In the Standard Model, i.e. in the absence of S- and T-type interactions, of = 1 and aar = —1/3. 
Any admixture of S to V (T to A) interaction in such a pure Fermi (Gamow- Teller) decay would 
result in a < 1 (a > —1/3). A measurement of a therefore yields information about the interactions 
involved. However, the neutrino cannot be detected directly and the 13—v angular correlation thus 
has to be inferred from other observables. From the properties of the general Hamiltonian of the 
weak interaction (Eq.Q and 0|) and of the Dirac 7-matrices it can be shown that the two leptons 
in /3-decay will be emitted preferably into the same direction for a V (T) interaction and into 
opposite directions for an S (A) interaction. This will lead to a relatively large energy of the recoil 
ion for a V (T) interaction and a relatively small recoil energy for an S (A) interaction (Fig.^. The 
WITCH experiment aims to measure the shape of the recoil energy spectrum with high precision 
to determine the (3 — v angular correlation parameter a and from this deduce a limit on a possible 
scalar or tensor admixture in weak interactions. 

Due to its feature of being able to measure the recoil ion energy spectrum in nuclear /?-decay, the 
WITCH experiment also provides interesting possibilities for other observables, e.g. one can also 
determine F/GT mixing ratios, Q- values, EC / (3'^ branching ratios and charge state distributions 
(1^ . [2^ . In addition, the ion cloud in the decay Penning trap can also be used for (3 and 7 
spectroscopy, this time not using the recoil spectrometer but with [3 and 7 detectors added to the 
WITCH set-up. For example, a (3 detector on the axis of the decay trap at some distance behind 
the trap opens the possibility for (3 spectroscopy with a pure sample without any scattering of the 
(3 particles in the source. Adding 7 detectors around the center electrode of the trap will make 
also 7 spectroscopy possible. 
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2 Experiment 



2.1 Principle 

An experiment to measure the recoil energy spectrum in nuclear /3-decay faces two major difficulties: 
(1) the /3-emitter is usually embedded in matter, which causes a distortion of the recoil ion spectrum 
due to energy losses caused by ion scattering in the source, and (2) the recoil ions have very low 
kinetic energy rendering a precise energy measurement difficult. 

In order to avoid the first problem and to be as independent as possible from the properties of 
the isotopes to be used, the WITCH experiment uses a double Penning trap structure to store 
radioactive ions. The ion cloud in the second trap, the decay trap, constitutes the source for the 
experiment, where the ions are kept for several half-lives, i.e. of the order of 1 to 10 s for the cases 
of interest 

To solve the second problem and measure the recoil energy spectrum a retardation spectrometer 
is used. The working principle of such a device is similar to the_5-spectrometers used for the 
determination of the neutrino rest-mass in Mainz 0| and Troitsk The WITCH spectrometer 
consists of two magnets: the first one providing a field B^ax = 9 T, the second one providing 
Bmin = 0.1 T, and an electrostatic retardation system. Recoil ions are created in the strong 
magnetic field region (i.e. in the Penning trap) and pass on their way to the detector the region 
with low magnetic field (i.e. the retardation section of the spectrometer). Provided that the fields 
change sufficiently slow along the path of the ions, their motion can be considered as adiabatic. 
According to the principle of adiabatic invariance of the magnetic fiux p\l B — const, where 
Pi_ is the momentum projection perpendicular to the magnetic field B. From this it follows that 
the radial kinetic energies of an ion in the trap {E*^°^^) and in the retardation section {E'^f^'^f) 
are related as El'll^J'l / E^j^,":^ = Brmn/Bmax- Thus a fraction 

1 - B„,,,jB„^ax ~ 98.9% (4) 

of the energy of the ion motion perpendicular to the magnetic field lines will be converted into 
energy of the ion motion along the magnetic field lines. The total kinetic energy of the recoil ions 
can be probed in the homogeneous region of low magnetic field Bmin by retarding them with a 
well-defined electrostatic potential. By counting how many ions pass the analysis plan e for different 
retardation voltages, the cumulative recoil energy spectrum can be measured 



2.2 Overview of the set-up 

The general scheme of the set-up can be seen in Fig. [2l The installation at the ISOLDE facility 
at CERN was recently completed. In a first step the ions produced by ISOLDE [11 El El get 
trapped and cooled by REXTRAP E3- soon as a sufficient amount of ions (viz. 10^ to 10^ 
ions) has been collected by REXTRAP they are ejected as a 60 keV (optionally 30 keV) bunch 
and are transmitted through the horizontal beamline (HBL) of WITCH into the vertical beamfine 
(VBL). There the ions are electrostatically decelerated from 60 keV to ~80 eV in several steps. 
In order to avoid a high voltage platform a pulsed cavity is used (l^ (named Pulsed Drift Tube, 
PDT, in the case of WITCH Q). When the ion bunch is inside the PDT the potential of the 
cavity is switched down over the range of 60 kV (30 kV) from 52 kV (26 kV) down to -8 kV (-4 kV). 
In this way the kinetic energy of the ions is not changed while the potential energy is shifted to 
-8 kV (-4 kV), so that the total energy becomes nearly zero (^80 eV in practice). The ions can 
then be captured in the cooler trap (the first Penning trap of the WITCH set-up), which is at 
ground potential. In this cooler trap the ion cloud is prepared (i.e. cooled and centered) before 
being ejected through the pumping diaphragm (which separates the vacua of the two traps) into 
the second Penning trap, the decay trap. The latter is placed at the entrance of the retardation 
spectrometer. After /?-decay the recoil ions emitted into the direction of the spectrometer spiral 
from the trap, which is in the strong magnetic field, into the weak field region. In the homogeneous 
low-field region the kinetic energy of the ions is then probed by the retardation potential (Fig. . 
The ions that pass this analysis region are re-accelerated to ^10 keV to get off the magnetic field 
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lines. The re-acceleration also ensures a constant detection efficiency for all recoil ion energies. 
Finally, the ions are focused with an Einzel lens onto the detection micro-channel plate (MCP) 
detector. For normalization purposes several /3-detectors are installed in the spectrometer section 
too (also to check the /?-simulations, see Sec. 12.611 The recoil spectrum can be measured 

setting one retardation step for one trap load (in this case a normalization, e.g. by counting the 
/^-particles, is necessary) or scanning all retardation steps during the same trap load (in which 
case a correction for the time-dependence of the count rate due to the isotope half-life has to be 
performed). The latter case also allows to avoid a possible effect of the MCP degradation on the 
shape of the recoil spectrum. 

2.3 Response function 

Since the WITCH experiment measures the recoil energy spectrum, good knowledge of the spec- 
trometer response function is of high importance. The response function has been investigated 
both analytically as well as with numerical calculations of the ion trajectories through the spec- 
trometer Two important issues regarding this response function are discussed below. 

2.3.1 Influence of the residual gas 

The residual pressure in the spectrometer section (and also in the decay trap) should be as low as 
possible. Indeed, in order to measure precisely the recoil spectrum it is very important to avoid 
scattering of the recoil ions. The residual gas can also cause another problem for the experiment, 
namely charge exchange which leads to neutralization of the ions such that they can not be probed 
anymore by the retardation principle. 

The influence of ion scattering was investigated via simulations of the ion trajectories, taking 
into account a Stokes force to describe the damping of the ion motion. The results are presented in 
Fig. 31 These simulations were performed for argon gas in both traps and in the spectrometer. The 
same pressure is used everywhere. For helium gas the indicated pressures have to be multiplied 
by a factor of « 8. One can see that with increasing rest gas pressure the response function 
broadens. It can also be noted that for a pressure better then 10~^ mbar the energy distribution 
deviates only a little bit from the ideal case p = mbar, meaning that for a successful recoil 
spectrum measurement the residual gas pressure has to be < 10~^ mbar in the decay trap and 
in the spectrometer. Up to now no estimates were made for the effect of charge exchange on the 
WITCH response function. 

2.3.2 Doppler broadening 

The ideal response function was derived assuming that the ions are at rest. The real situation, 
however, differs since the ion cloud in the decay trap is at least at room temperature. The velocity 
of the recoil ion is thus the superposition of the velocity of the mother nucleus in the trap and the 
velocity obtained due to the beta decay. This definitely affects the recoil spectrum. Assuming a 
MaxweUian ion velocity distribution and taking into account that the energy at room temperature 
(0.025 eV) is much smaller than the recoil energy 0(100 eV), it can be shown ^3 that the energy 
shift e due to Doppler broadening results in the Gaussian distribution: 



where E^f^ is the recoil ion energy obtained in the beta decay, k is the Boltzmann constant and T 
is the ion cloud temperature. Eq.© yields, for e.g. E^'r^"'' = 280 eV and T = 300 K, cr, = 3.8 eV. 
This means that the mono-energetic response function broadens significantly if the ion cloud is at 
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room temperature. To get the real response function for the WITCH spectrometer one has to fold 
the ideal mono-energetic response function with this broadened energy distribution. Fig. [3 shows 
the WITCH response function thus obtained for different ion cloud temperatures. As one can see, 
down to liquid nitrogen temperature (i.e. 77 K) the shape of the response function is defined by 
Doppler broadening and only at the temperature of liquid helium (i.e. 4 K) it approaches the 
shape of the ideal response function. It would thus be interesting to consider to cool (at least to 
77 K) both Penning traps. 

The response function can be determined experimentally by measuring recoil ions from electron 
capture decay (EC), which leads to a mono-energetic peak that is, if it accompanies /?"'" decay, above 
the endpoint energy of the continuous recoil spectrum from /3+ decay (Fig. see Sec. I2.4|l . By 
measuring these EC peaks for several suitable nuclides with different decay energies, an energy 
calibration of the spectrometer can be performed as well. 



2.4 Recoil spectrum 

Electron shake-off after the /?-decay j^^l will cause the daughter ions to have different charge states 
q = n-e. The recoil ions with energy -E^f™*' and charge q will appear in the measured spectrum at a 
retardation voltage Uret = E^tn^^ ll due to the retardation principle. The measured spectrum will 
thus be a superposition of the spectra of the various charge states, each with different endpoints 
Uon = ^^"'^P°™*/ (n-e), n > 1, where jg ^^e endpoint energy of the recoil energy spectrum 

for n = 1 (Fig.EJ. Consequently, when measuring the full recoil spectrum up to the endpoint energy 
j^endpoint ^ the Upper half of the spectrum will consist purely of events from charge state n = 1. 
This upper half of the recoil spectrum is therefore the most interesting part for analyzing the 
spectrum shape • Note that the EC peaks of different charge states n for the same nuclide will 
also appear in the spectrum at different energies EEc/n, n>l. 

It is also planned to check the dependence of the shake-off probability on the recoil ion energy, 
as has been seen in ^He /3~ decay (l^ . This dependence is expected to be even larger in /3+ 
decay jl^l- This effect to first order distorts the recoil ion spectrum by (1 + s • Erecou)- The idea is 
to use the upper half of the l"*" charge state spectrum to fit simultaneously s and (3 — ly correlation 
coefficient a. A possible dependence of the shake-off process on the recoil energy can be checked 
by fitting spectra obtained for different charge states to a (not including the s parameter). 



2.5 Achievable precision 

A series of random integral recoil spectra have been generated from which the (3 — v angular 
correlation parameter a was fitted. The response function of the spectrometer was approximated 
by Gaussian with a —1%. No /3 background and no dependence of the shake-off on the recoil 
ion energy were considered. By changing the fit interval and the bin width the infiuence on the 
achievable precision for a was studied. These simulations (Fig.O show that, to reach a precision of 
Aa = 0.005, the total number of events in the differential energy spectrum should be iV = 10''-^ 10^ 
and a minimum of ng = 20 channels (i.e. retardation steps) in the upper half of the spectrum seems 
to be sufficient. Taking iV = 10^, tt-q = 20, the number of ions in one trap load Nioad = 10^ and 
the efficiency parameters for a fully optimized set-up (Table QJ, one can find that the measurement 
time needed to reach a precision of Aa = 0.005 (Fig. is tmeas ~ 3.6 days for "^^Ar (ii/2 = 1.78 s) 
and t„ieas ~ 13 days for ^Sm^j {ti/2 ~ 6.35 s) In fact, the results of the simulations presented 
in Fig. 13 show that the required precision of Aa = 0.005 might be achieved already with TV = 10^ 
which means that the measurement time reduces correspondingly, i.e. tmeas ~ 8.6 hours for '^^Ar 
and t„ieas ~ 1-3 days for ^^^Al. 



2.6 GEANT4 /^-particle simulations 

As was already mentioned in Sec. 12.21 normalization detectors (to count /3-particles) are available 
in the set-up too. Simulations with the GEANT4 simulation package |25l |2a| were performed to 
find the best suited place for installing these detectors in order to provide good statistics for the 
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normalization. The scheme of the simulated set-up is presented in Fig.|Hl(iop). The double Penning 
trap structure is described as a set of two solid copper cylinders (inner diameter 40 mm) separated 
by a copper pumping diaphragm with an opening of 4 mm diameter. All spectrometer electrodes 
as well the walls of the lower and upper magnet bore tubes are made from stainless steel. For the 
MCP detector a very simplified model is used, i.e. a solid disk of quartz material (Si02), 80 mm 
in diameter and 1 mm thick. 

The program was run for the case of '^^Ar, with an isotropic /3-particle emission distribution 
and a 5 mm radius ion cloud in the center of the decay trap. The positions and elements of 
the set-up which are hit first by the /3-particles after they leave the decay trap were determined. 
This is shown in Fig. |S| One can recognize the shape of certain electrodes and also the pumping 
diaphragm (at nearly -400 mm). The largest number of hits is in the pumping barrier: nearly 80% 
of the betas going in this direction hit the barrier. The second most hit element is the SPDRIFOl 
electrode. Taking into account that the size of the pumping barrier is much smaller than the size 
of the SPDRIFOl electrode, one can conclude that the pumping barrier is the most suitable place 
to install the /^-detectors. 

Another important result of these simulations is the estimate for the /3-background on the 
recoil ion detector. This MCP detector is introduced in the program as a very primitive model. 
Nevertheless, this allows to check how many particles will arrive on the detector. It appears that 
this amount is actually rather large: from a total of 250 000 simulated events 8608 betas arrive on 
the 8 cm diameter MCP detector. This number has to be compared to the number of ions arriving 
on the same detector. Considering the efficiencies for an ideal set-up (TableOJ it can be estimated 
that about 10 000 ions will arrive on the MCP detector. However, it might be that an MCP of 
4 cm diameter is already sufficient for detecting all ions. Since the /3-background events equally 
cover the surface of the MCP, a two times smaller diameter of the detector means a reduction of 
the /9-background by a factor of 4, corresponding to 2152 beta-particle events. This estimate of 
the /3-background is not yet complete however, as one still has to take into account the different 
detection efficiency of the MCP for ions of 0(10 keV) energy (the recoil ions are accelerated onto 
the MCP, see Sec. [221 and /3-particles of O(MeV) energy. The model used for the MCP detector 
is too primitive to answer this question. Measurements performed by the TRIUMF and Berkeley 
groups l^^l^ l showed that this efficiency is actually very close to the one for ions (i.e. the absolute 
MCP detection efficiency is about 50-f-75% for /3-particles of O(MeV) energy, compared to about 
60% for ions) . This means that if one assumes the MCP registration efficiency for betas and ions 
to be both 60%, the MCP detector will register ^1290 /3— particles against 6 000 ion events in the 
full recoil ion spectrum. 

The first estimates [i^l show that in order to achieve Aa = 0.005 with such (3 background an 
increase of the measurement time by about one order of magnitude might be necessary. However, if 
one considers that lO'' counts in the differential spectrum (not including (i background) is sufficient 
(Fig-El Sec. l2.5|l . the required precision on the P—v correlation (Aa = 0.005) might still be achieved 
in a realistic time period. These simulations, for sure, have to be checked experimentally and the 
issue of the MCP detection efficiency for ions and betas still has to be investigated in more detail 
too. 

3 Tests performed 

The WITCH set-up was completed only recently and commissioning was started already during 
the last phase of the installation. In most of the tests the off-line REXTRAP ion source (delivering 
a ■^^K ion beam) was used. The tests performed are discussed in this section. 

3.1 HBL 

The function of the HBL is to transfer the ion beam from the REXTRAP set-up into the WITCH 
vertical beamline as efficient as possible. Since the ISOLDE facility (and respectively, REXTRAP) 
is operated usually at either 60 kV or 30 kV, tests were done using both high voltage settings. The 
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tests performed showed that the tuning can be done in such a way that there is no significant loss 
of beam intensity through the HBL, i.e. nearly 100% transmission efficiency is obtained. To tune 
the voltages and monitor the beam, the three diagnostics in the HBL and the first one in the VBL 
are used. 

3.2 VBL 

The main purpose of the VBL tests that were performed till now was to prove the functionality of 
the PDT, check its efficiency and investigate and optimize the injection of the ion beam into the 
WITCH magnetic field. 

From the values for the resistances and capacitances in the HV switch box scheme (Fig. 
the HV switching time constant tpdt was estimated to be 0.18/-ts. This provides a pulse down 
time tswitch of the PDT of about 1.3 //s (1.2 //s) for a 60 kV (30 kV) ISOLDE beam, with tswitch 
being defined such that the energy of the ions after the switching differs less than 50 eV from the 
required value The travel time of 30 keV ^^K"*" ions in the PDT (for a combination of HV = 
+26/-4 kV) is ~ 5^s. This gives an ~ 3.7/is time window to pulse down the ion beam. However, 
the ion bunch ejected from the REXTRAP set-up has a time structure that is longer than this time 
window (the typical bunch length is ^ 10 ^s), meaning that it is not possible to pulse the complete 
bunch but at best only about 43% of it for a combination of HV = +26/-4 kV [l^l- The original 
bunch can thus be considered to consist of different parts: 1) ions which will pass through the 
PDT before the pulsing down starts, 2) an intermediate part containing partially bunched ions, 3) 
well bunched ions, 4) another intermediate part of partially pulsed down ions and, finally, 5) ions 
which enter the PDT when this is already at low voltage (Fig. top left corner). For the cases 
1) and 5) these ions will have ~ 60 keV (30 keV) after the PDT, i.e. they are much faster than 
the well-bunched ions. This leads to the time structure of the signal after the PDT shown also 
in Fig. EH In this figure the simulated (using the SIMION simulation routine j2^) and measured 
spectra are compared. As can be seen, good agreement between simulations and the measurement 
is obtained. 

The efficiency of the PDT, 7]pdt, was measured using two diagnostic MCP detectors: one in 
front of the PDT and one behind it at ^^40 cm from the exit of the PDT. The resulting efficiency 
Vpdt = 8% is less than the expected value of ~43% but is of the right order of magnitude 

From the signal of the MCP detector behind the trap structure (i.e. in the magnetic field) 
the overall efficiency of the vertical beamline, including the efficiency of the pulsed drift tube, 
was found to be between 0.1% and 1% (for these measurements ions were not trapped but only 
sent through the traps) (Table . Additional measurements were performed in the mean time 
with improved MCP diagnostic detectors j^^, i.e. with several Ni meshes in front to reduce the 
incoming beam intensity, and a grid-anode to determine the beam size. These have revealed a 
significant decrease of beam quality after the PDT. This is being looked into in more detail now. 

3.3 Traps 

In the initial phase, i.e. for a simple mode WITCH operation and tests of the spectrometer, a 
sophisticated use of the traps is not required. It was therefore not our aim to perform already now 
a detailed investigation of their properties. Rather we wanted to verify the simple operation of 
the traps, i.e. try out basic trap mechanisms, understand the behavior of ions in the traps and, if 
possible, optimize the parameters in order to find a suitable trapping regime. 

3.3.1 Buffer gas cooling 

An important step in verifying the operation of Penning traps is to check the buffer gas cooHng 
procedure since this is usually a first main step in the ion cooHng. The buffer gas used is high 
purity ^He (quality: Helium 57 or >99.9997%), but the transfer line is not yet equipped with a 
cold trap or other purification system. The buffer gas pressure measurement is done at a point 
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after the gas dosing valve (Pfeiffer Vacuum, RME005) and before the transfer line enters into the 
vacuum chamber. With this installation the buffer gas pressure is regulated via a feedback loop. 

A scan of the MCP signal as a function of the cooling time was performed for a buffer gas 
pressure of 5 mbar at the gas dosing valve position. "^^K"'" ions are trapped in the cooler trap, 
cooled there for some time and then sent (without capturing) through the decay trap onto the 
MCP detector behind the WITCH trap structure. It was found that after ~200 ms of cooHng 
the MCP signal splits in two different peaks (Fis. This means that other ions then '^^K"'" are 
present too, for instance because the buffer gas is not clean enough. Using dipole excitation of 
the reduced cyclotron motion {v+), the first peak was identified as mass 19 (H2OH"'') while the 
second one corresponds to '^^K"'". Qualitatively the effect of the cooHng can be seen in Fig. ITTl 
after 100 ms of cooHng the TOF position of the '^^K+ peak is ~ 70 while after 200 ms of cooHng 
the peak appears at ~ 100 fis. 

3.3.2 Excitations 

In Fig. El the oscilloscope pictures corresponding to the dipole iy+ excitation are shown. With 
no applied excitation two peaks (H2OH+ and '^^K+) are visible. Next, the foHowing scheme is 
used: first the ions are cooled by colHsions with the buffer gas atoms in the cooler trap during 
200 ms, thereafter they are excited at v+{^^K) = 2367100 Hz (in a 6 T B-field) for 100 ms with 
an amplitude of A^^ = 2 V and, finally, they are extracted fFig. \12b). As can be seen the '^^K"'" 
peak disappears, i.e. with v+{^^K) excitation the '^^K"'' ions were brought to a radius larger than 
the radius of the pumping diaphragm. In another test the ions in the first peak were excited at 
the reduced cyclotron frequency i^rf = 4845000 Hz ~ z^+(mass 19) {Ai,^ = 0.5 V) during 100 ms 
(in this case the ions were first cooled during 400 ms) (Fig. Efc)- The same effect is obtained, 
i.e. the corresponding peak disappears. The increase of the MCP signal after removing the first 
peak could be related to MCP effects (see Sec. K-i.5|l . No systematic scan of the number of ions 
ejected from the trap as a function of the excitation frequency was performed as yet, but only a 
qualitative study of the TOF osciHoscope spectrum. These tests nevertheless show that the 
excitation works for WITCH and aHows to separate different masses. 

A mass selective removal of unwanted species can be achieved via a combination of dipole 
excitation of the magnetron motion (i^_) and quadrupole excitation at the true cyclotron frequency 
(vc) This technique was tried qualitatively for '^^K+ ions. A dipole ly- excitation (at v- = 

130 Hz, determined in a similar way as v+, i.e. checking at which frequency all ion peaks disappear 
from the osciHoscope spectrum) with amplitude Ai,_ = 150 mV was applied for 50 ms, and followed 
by a quadrupole Vc excitation (i^c('^^K) = 3553729 Hz) with amplitude A^^{^^K) — 1.6 V for 3 ms 
(B-field is 9 T). As in the previous tests no systematic scan of the excitation frequencies was 
performed but only a visual analysis of the TOF oscilloscope spectrum. The corresponding steps 
of the process are shown in Fig. El Fig. [T3I 7, displays the situation before any excitation. The 
second peak corresponds to •^^K+. When a dipole excitation is applied, all ions are driven out 
as can be seen in Fig. \Wh (both peaks disappear). If now a quadrupole RF-field at frequency 
i^c('^^K) is used, one can expect the '^^K^ ions to recenter while other impurities should disappear. 
However, while it is clear from Fig.ll3fc that there is indeed no other species present than '^^K"'", the 
signal corresponding to the '^^K+ ions is significantly smaHer and broader than the one without any 
excitation. A possible reason for this can be either wrongly chosen parameters (so that the ions hit 
the electrode and are lost) or some electronics problem (e.g. electronic noise in some channels, a 
different capacitance of the track: connection wires + electrodes, or imperfection of power supplies 
used). Also, the behaviour of the ions might have been infiuenced by space charge effects. 

The dipole v- excitation was also applied while working with ^^Ne"*" ions and the frequency 
used there was ~ 140 Hz (to successfuHy remove aH species; B-field was 9 T). Combining this 
value with the one in the '^^K test leads to = 135(5) Hz. Using now this value of and an 
estimation of the magnetic field in the cooler trap center (see next Sec. 13. 411 . one can deduce the t rap 
characteristic parameter Uo/(P = 1.53(6) • 10^ V/m^ for WITCH. For the ISOLTRAP [alE3Jli 
cooler trap, which is very similar to the WITCH cooler trap, one has Uo/(P = 1.8- 10^ V/m^ |34j . 
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3.4 Magnetic field 



The necessity to know precisely the magnetic field at the trap center is based on the following two 
factors: 1) the cyclotron frequency, and therefore the centering and cooling of the ion cloud as 
well as the mass-selectivity, are directly determined by the value of the field, while 2) the response 
function of the WITCH spectrometer also depends on it. Originally the magnetic field map was 
provided by Oxford Instruments for both magnets separately. However, this field map resulted 
from a calculation/approximation based on measurements that were made at the factory, prior to 
the delivery and installation of the system at CERN. Also, inserting the traps structure with all 
the wiring may change the field strength in the trap centers due to the magnetic susceptibility of 
the materials used. Direct measurement of the field with an NMR probe is hardly possible because 
of the very difficult access to the area. 

An elegant way to estimate the magnetic field is to use a known isotope (i.e. with known 
mass m) and experimentally find the proper cyclotron frequency Vc. This yields enough infor- 
mation to determine the magnetic field B From the excitation tests performed with a 
^^K+ beam, j/c('^^K) was found via quadrupole excitation of ^^K+ (direct determination, Vc = 
3553729(2000) Hz) but also via dipole excitation oi^^K+ (which gives only v+ = 3554019(2000) Hz; 
in another test = 135(5) Hz was found, see plHJ. The estimated uncertainty for and v+ is 
based on two measurements of v+{^^K) at 6 T field: v^^\^^K) = 2363170 Hz and v'^^\^^K) = 
2367100 Hz (FigEll, which leads to Av+{^^K) ^ 2000 Hz. The final result for vJ(^^K) is 
Dc = 3553900(1400) Hz, which corresponds to a magnetic field in the center of the cooler trap 
Bcooier = 9.018(4) T (according to the field map from Oxford Instruments the field in the center 
of the cooler trap is 9.0011 T). The set field was 9 T for the lower magnet and 0.1 T for the top 
magnet of the system (the WITCH magnet system allows to set any possible combination of the 
magnetic fields in the range < Bmax < 9 T and < Bmin < 0.2 T). 

3.5 MCP regime 

The working principle of an MCP detector can be found in e.g. Under certain conditions 

an MCP is not sensitive anymore to the number of incident particles (Fig.^Jl, meaning that the 
detector misses some events. This state is known as the saturation of the MCP. The typical dead- 
time of one MCP channel is in the order of several tens of milliseconds j^. For the WITCH 
diagnostic MCPs it is ^^30 ms. The saturation of an MCP depends on both the ion current density 
(i.e. the number of ions per MCP channel and per second) and the MCP acceleration voltage. The 
effect of the ion current density is shown in FigEI when the number of ions in the bunch exceeds a 
certain value the MCP signal remains constant. However, already much earlier the dependence of 
the MCP signal on the number of incident ions ceases to be linear. This behavior can be explained 
by partial saturation: Fig. IT^ shows that at MCP HV=1.25 kV the ^^Na"*" signal has a block 
shape while at higher MCP voltages the total signal still increases but shows a significant drop in 
intensity for later arriving ions. This means that the early arriving ions saturate a certain fraction 
of the MCP channels, leading to a decrease of the MCP registration efficiency for later ions. This 
can also be seen from Fig. El when the early arriving ions are removed with a time window before 
they reach the detector (the WITCH beam gate is used for this), the signal corresponding to the 
late ions increases. 

The effect described above infiuences the measurements and has to be taken into account for 
efficiency estimates and during the beam tuning. Saturation of the diagnostic MCP's can be 
avoided by adding several Ni meshes in front and using appropriate MCP acceleration voltages. 
The dead-time and saturation of the MCP can also infiuence the recoil spectrum measurement, 
since in real measurement conditions the recoil ions are supposed to reach the recoil MCP detector 
at a rate > 10^ Hz. Results of a careful study of the MCP response to high intensity pulsed beams 
are described in |3a|. 
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3.6 First radioactive ions 



In November 2004 WITCH got its first radioactive beam time (with 35Ar+). A CaO ISOLDE 
target and a plasma ion source with cold transfer line were used. The rate for this particular 
run was about 5 • 10^ atoms/s, being somewhat lower then expected. Coupled with the fact that 
the efficiency of WITCH is not yet high enough to deal with this intensity (see Tabled a recoil 
energy spectrum measurement was not yet possible. Increase of ISOLDE rates of 40 times can be 
expected based on past demonstrated yields For WITCH ongoing optimization should yield 
an increase in overall efiiciency of a factor of ten or more. Still, the decay of '^^Ar was observed 
on the first VBL MCP detector (Fig. El. The half life of ^^Ar obtained as the weighted average 
of two short measurements is: Ti/2('^^Ar) = 1.70(5) s (the value in literature is 1.775(4) s). This 
showed that ISOLDE delivered a clean '^^Ar+ beam with no radioactive contaminant, as is required 
for the planned recoil spectrum measurements. 

4 Outlook & improvements 

The WITCH set-up was completed and first commissioning tests performed only recently. There is 
still room for many improvements and more tests are necessary to better understand the behaviour 
of the different parts of the set-up. 

The non-pulsed high energy 60 keV (30 keV) ions arrive first on the diagnostic MCP and can 
cause saturation of the detector, reducing its sensitivity, i.e. disturbing the tuning and efficiency 
measurements. Part of these energetic ions also reaches the detection MCP at the end of the 
spectrometer, in spite of the magnetic field. This will infiuence the measurement of the recoil ion 
spectrum. Another drawback is that during a radioactive run, the decays of non-pulsed 60 keV 
(30 keV) radioactive ions implanted directly on the detection MCP will lead to additional back- 
ground. All these problems can be avoided if one uses a beam gate installed in the HBL (or VBL) 
in order to select only the part of the original beam corresponding to the correctly pulsed ions. 
The required electronics to switch the voltages in the range of 1000 V within several 100 ns is 
currently being developed. 

A new system of VBL diagnostics is currently being prepared. It is based on an MCP detector 
with split anode and a Ni-mesh in front of it. The latter reduces the intensity of the incoming 
beam in order to avoid saturation of the detector. The transparency of this mesh can be measured 
to good precision with laser light. The split anode system provides the possibility to check the 
beam size and its position "on-line". The combination of the Ni-mesh and the split anode will allow 
to avoid problems caused by the saturation of the MCP during the beam tuning. 

Based on the new HV switch system of ISOLTRAP a new HV scheme was developed in close 
cooperation with a company^. This new system is more reliable for 60 kV switching and has in 
addition the advantage that the switching time is improved by a factor of 2 3. This is done with 
an approach of a clamping diode, which ties the decreasing voltage to a pre-defined value for a 
limited period of time. The switching process in this case starts as usual ('Fig. ll8|l . i.e. the voltage 
of the PDT goes down towards the negative biasing voltage. After roughly 600 ns, corresponding 
to 3 X TpDT, the voltage has dropped below the value of an auxiliary voltage supply to which the 
diode is connected. The diode therefore becomes conducting and prevents the PDT-voltage from 
a further decrease. 

The trap tests showed that the buffer gas of the cooler trap contains impurities. The main 
effect of this is that the ions of interest can be neutralized via charge exchange and can thus 
escape from the trap. To avoid this, the external gas fine has to be as short as possible (the 
internal part cannot be changed), while in addition it has to be cleaned, baked and pumped to 
remove contaminations. To further clean the buffer gas one can in addition install a cold trap or 
use a commercially available He purifier. During the commissioning period it was realized that a 
higher energy for the ions leaving the PDT improves the injection in the cooler trap. However, 
this requires an upgrade of the end cap power suppHes. The corresponding electronics is currently 

^ Dr. Stefan Stahl - Elektronik-Beratung, Sonderanfertigungen ■ Kellerweg 2.3, D - 67.582 Mettenheim ■ Germany. 
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being developed. With respect to the detection part, an 8 cm diameter MCP detector with position 
sensitive anode will be used to study the size of the recoil ion beam, the possible dependence of the 
beam size on the ion energy, as well as the /j-background. This MCP and the necessary electronics 
will be provided by the LPC-Caen group |4lll42l|. 

To carry out efficiency tests and improve the beam tuning in WITCH, it is necessary to transport 
ion beams through the complete WITCH beamhne. The REXTRAP ion source is often needed by 
the REXTRAP team for tests and, in addition, can not be used during experiments involving the 
REXTRAP set-up (since this ion source blocks the ISOLDE beam at the entrance of REXTRAP). 
A design study was therefore started to develop an ion source for WITCH similar to the one of 
REXTRAP and to implement this in the horizontal beamline. 

The vacuum of the WITCH system is at present reasonably good 10^^ -f- 10"'' mbar) for 
normal WITCH operation but can still be improved to avoid pressure related systematic effects 
and to reduce the charge exchange probability. The WITCH spectrometer was designed with the 
possibility to use non-evaporable getters (NEG) and it is planned to put this system in operation 
soon. 

Finally, additional tests are being prepared in order to improve the WITCH efficiency. For 
instance, one may try to optimize REXTRAP operation so as to reduce the ion bunch length, try 
to determine the rest gas pressure in the WITCH traps, study the size of the ion cloud, investigate 
space charge effects, and try to obtain an optimal cooling time. In addition, investigation of the 
spectrometer and its response function, as well as measurements to check the /?-background on the 
main detector are planned as well. 

5 Conclusion 

The installation period of the WITCH set-up, which was developed over the last few years, was 
finished in autumn 2004 while intensive commissioning of the set-up was performed during the 
whole year 2004. The main aim of these tests was to check the operation of the beam transport, the 
pulsing down of the ion beam and the injection of ions into the high magnetic field, to test the trap 
basics and check the spectrometer operation, and finally, to optimize as many settings as possible. 
These tests showed that the full set-up up to the spectrometer is now operational, although several 
efficiencies still have to be improved. The present overall efficiency of the experiment did not allow 
to actively test the retardation spectrometer which can only be done fully with radioactive ions. 
The results of the commissioning stage were carefully analyzed and possible improvements were 
suggested. This includes both technical modifications as well as the necessary tests in order to 
optimize the set-up and achieve the required efficiency. 
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Figure 2: Schematic view of the WITCH set-up. The abbreviations used are: HBL horizontal 
beamline, VBL vertical beamline, PDT pulsed drift tube. The HBL is 90° rotated (i.e. the top 
view of the HBL is shown). 
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z (cm) 

Figure 3: Magnetic field (solid line) and retardation electric field (dashed line) profile on the axis 
of the spectrometer (z-axis). The retardation field has been calculated for Uret = 100 V. z=0 
corresponds to the center of the 9 T magnet. The positions of the traps, the analysis plane and 
the recoil ion detector are also indicated. 
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Figure 4: Influence of difi^erent rest gas pressures on the response function of the WITCH spectrom- 
eter: mbar (soUd Une), 10~^ mbar (dashed Une), 10~^ mbar (dash dotted hne) and 10~^ mbar 
(dotted hne). The simulations are performed for a particle with ErecoU = 500 eV and mass=50 amu 
in argon gas. 
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Figure 5: Doppler broadening of the response function of the WITCH spectrometer for different 
temperatures of the ion cloud. Calculations performed for ErecoU = 280 eV without taking into 
account the rest gas pressure. 
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Figure 6: Differential recoil spectrum calculated for the /3+-decay of "^^"K. The charge distribution 
for the /S^-decay is taken from [iJl, while the charge state distribution for electron capture is 
unknown for this isotope. In order to still get a qualitative idea the known charge state distribution 
for '^^Ar is therefore shown (from 
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Events in diff. spectrum 

Figure 7: Estimated precision on the 13— v angular correlation coefficient a (for ^^™A1, C.L.=68.3%) 
as a function of the total number of events N in the differential recoil spectrum when three different 
energy intervals near the endpoint (at 280.6 eV for ^^™A1) are considered for analysis. 
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Figure 8: Plot of the points where /^-particles hit the set-up for the first time. '^^Ar spectrum, 5 mm 
cloud, isotropic distribution. The center of the decay trap is at -323 mm. On top the WITCH 
set-up in the GEANT4 simulation program is shown. The different spectrometer (SP) electrodes 
are indicated. 
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Figure 9: Electrical scheme of the HV switch system for 60 kV. 
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Figure 10: TOF spectra of ^^K after the PDT. The simulation is for a HV switch time constant = 
0.2 ^s. The measured spectrum is the MCP signal inverted and scaled to the simulated spectrum. 
Measurement and simulation are both for HV PDT = 26 kV. The zero of the TOF -axis corresponds 
to the start of the HV switching. In top left corner the schematic of the original ion bunch is shown: 
1) "early" ions, 2) partially bunched "early", 3) well bunched ions, 4) partially bunched "late", 5) 
"late" ions, a, b, c and d correspond to the transitions between different parts. 
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Figure 11: Oscilloscope TOF pictures observed with an MCP for different cooling times (100 ms, 
200 ms, 500 ms) in the cooler trap. The buffer gas pressure was 5 mbar (at the gas dosing valve 
position). 
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Figure 12: Identification of different peaks from the TOF observed with an MCP detector after 
different excitations in the cooler trap. Ions are a) first cooled for 200 ms, without exciting them 
and then b) dipole excited at Vrf ~ v^i^'^K) = 2367100 Hz {A„^ = 2 V) during 100 ms and c) 
excited at the reduced cyclotron frequency Vrf = 4845000 Hz ~ i^+(mass 19) {A^^ = 0.5 V) during 
100 ms (in this last case the ions were first cooled during 400 ms). The peak corresponding to 
the excited mass disappears. The magnetic field was 6 T. The He buffer-gas pressure was 5 mbar. 
The increase of the MCP signal after removing the first peak (case (c)) could be related to MCP 
effects (see Sec. I3.R|I 
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Figure 13: Mass selective cooling of '^^K: a) no excitation is applied 6) excitation; all peaks 
disappear c) quadrupole i^c('^^K) excitation; '^^K+ ions come back but the signal amplitude is, for 
some as yet unknown reason, much lower. 
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Figure 14: Summed MCP signal as a function of the number of incident particles (per bunch) for 
MCP HV = 1.4 kV (third VBL diagnostics, the MCP sensitive area is 18 mm diameter). The plot 
shows a combination of two independent measurements. As can be seen, the MCP saturates if the 
number of incident particles is > 2 • 10^ per bunch. 
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Figure 15: Change of MCP signal shape with applied HV, and the saturation effect (measured 
with third VBL diagnostics). 
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HV=1.6 kV): removing early arriving ions with the beam gate increases the MCP signal of later 
arriving ions. 
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Figure 17: ''^Ar half-life measurement on the first VBL diagnostic MCP detector. 
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Figure 18: HV of the PDT as a function of time in case of a standard exponential decrease with 
TpDT = 0.2 fj.s (1) and using the clamping diode (2) (from J4(\). 
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Description oi the parameter 


Efficiency 




ideal set-up 


achieved 


HBL efficiency, 'qHBL 


100% 


-100% 


PDT efficiency, tjpdt 


43% 


8% 


Injection into the magnetic field, ^injection 


100% 


1%-10% 


Trapping in the cooler trap 


100% 


-60% 


Losses during cooling 


100% 


_75%a),b),c) 


Efficiency of transfer between traps 


100% 


-80% 


Losses in the decay trap 


100% 


100% "'^'"^ 


Fraction of ions leaving the decay trap, taking 


40% 


not yet studied 


into account the cut-off angle (eg) and solid angle 












Shake-off for charge state n— I, P(n=i) 


10% 


not yet studied 


Transmission through spectrometer 


100% ''^ 


not yet studied 


MCP efficiency, cmcp 


60% 


52.3(3)% 


Total efficiency 


^ 1% 


-0.6- (10-3^10-2)%/) 



"'' 100% means that there are no losses. 

after 300 ms of cooling, with 5 mbar He buffer gas pressure at the measurement position. 

estimated from the tests performed, but one of the other effects (e.g. MCP regime) was not taken into 
account. 

ion energy dependent value. 

from E. Lienard et al., Nucl. Instr. and Meth. A, 551 (2005) 375. 

when an efficiency has not been determined yet the values for an ideal set-up are taken. 
Table 1: Efficiencies for a fully optimized WITCH set-up and the currently achieved values. 
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